Meso-scale turbulence in living fluids 
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Turbulence is ubiquitous, from oceanic currents to small-scale bi- 
ological and quantum systems. Self-sustained turbulent motion in 
microbial suspensions presents an intriguing example of collective 
dynamical behavior amongst the simplest forms of life, and is impor- 
tant for fluid mixing and molecular transport on the microscale. The 
mathematical characterization of turbulence phenomena in active 
non-equilibrium fluids proves even more difficult than for conven- 
tional liquids or gases. It is not known which features of turbulent 
phases in living matter are universal or system-specific, or which 
generalizations of the Navier-Stokes equations are able to describe 
them adequately. Here, we combine experiments, particle simula- 
tions, and continuum theory to identify the statistical properties of 
self-sustained meso-scale turbulence in active systems. To study 
how dimensionality and boundary conditions affect collective bacte- 
rial dynamics, we measured energy spectra and structure functions 
in dense Bacillus subtilis suspensions in quasi-2D and 3D geome- 
tries. Our experimental results for the bacterial flow statistics agree 
well with predictions from a minimal model for self-propelled rods, 
suggesting that at high concentrations the collective motion of the 
bacteria is dominated by short-range interactions. To provide a ba- 
sis for future theoretical studies, we propose a minimal continuum 
model for incompressible bacterial flow. A detailed numerical analy- 
sis of the 2D case shows that this theory can reproduce many of the 
experimentally observed features of self-sustained active turbulence. 
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Abbreviations: SPR, self-propelled rod; Re, Reynolds number 

Simple forms of life, like amoebae or bacteria, self-organize 
into remarkable macroscopic patterns [TJ [2] , ranging from 
extended networks [3 [4] to complex vortices [5l 161171 151 l9l [TO] 
and swarms [TT]. These structures often bear a striking re- 
semblance to assemblies of higher organisms (e.g., flocks of 
birds 12 or schools of fish [131 114]), and present important 
biological model systems to study non-equilibrium phases and 
their transitions 15, 16, 17 . A particularly interesting man- 
ifestation of collective behavior in microbial suspensions is 
the emergence of meso-scale turbulent motion 7, 8, 18l 119] . 
Driven by the microorganisms' self-propulsion and their mu- 
tual interactions, such self-sustained 'active turbulence' can 
have profound effects on nutrient mixing and molecular trans- 
port in microbiological systems 2, 20, 21, 22 . However, in 
spite of recent progress [19] [23] [24] [25], the phenomenol- 
ogy of turbulent bacterial dynamics is scarcely understood, 
and a commonly accepted theoretical description is lack- 
ing [2] 1161 126] . The latter fact may not be surprising given 
that a comprehensive mathematical characterization of tur- 
bulence in conventional fluids has remained elusive after more 
than a century of intense research [27] . 

In view of the various physical and chemical pathways 
through which bacteria may communicate [T] 1111 128] , a basic 
yet unsolved problem is to identify those interactions that are 
responsible for the emergence of collective behavior in dense 
suspensions [2j 1291 [30] . Answering this question is essential 
for understanding whether physical mechanisms such as flag- 
ellar bundling or hydrodynamic long-range interactions are 



relevant to collective bacterial motion; it is also crucial for 
constraining the vast number of theoretical models that have 
been proposed during the past two decades [2] [16] [19] [31] [32] , 
but have yet to be tested against experiments. An equally 
important, unresolved issue pertains to the 'universality' of 
turbulent phenomena in active systems and their relation to 
turbulence in passive fluids [27]. In ordinary liquids and gases, 
such as water or air, turbulent vortices form due to external 
forcing if the Reynolds number, the ratio of inertial to viscous 
forces, is very large (Re ^> 1). By contrast, bacteria provide 
an internal microscopic forcing and operate at Re ~ 10 -5 [33] . 
It is therefore unclear how, or to what extent, the characteris- 
tics of self-sustained turbulent states in microbial suspensions 
differ from those of classical turbulence in passive fluids. 

Here, we combine numerical simulations, high-speed mi- 
croscopic imaging and continuum theory to identify generic 
statistical properties of active turbulent motion in dense bac- 
terial systems, using Bacillus subtilis as a model organism. 
Unlike previous investigations of collective bacterial swimming 
in 2D free-standing films [8] and 3D bulk suspensions with 
liquid-gas interfaces [7] [2H [25], we conducted experiments 
in closed quasi-2D and 3D microfluidic chambers to mini- 
mize external influences and to compare the effects of bound- 
ary conditions and dimensionality. Our analysis focusses on 
traditional turbulence measures, such as energy spectra and 
velocity structure functions [27] [34] [35]. These quantities 
have been widely studied for turbulent high- Re Navier-Stokes 
flow [27] [36] [37] [38] [39] [40] [41], but their characteristics are 
largely unknown for active fluids. We compare our experi- 
mental results with large-scale simulations of a 2D minimal 
model for self-propelled rods (SPRs). In the past, similar 
models [42] have proven useful for identifying generic aspects 
of flocking and swarming in active systems 43 , 44 . We find 
that, although the SPR model neglects details of bacterial 
cell-cell interactions, it is able to reproduce many features of 
our experimental data, suggesting that collective bacterial dy- 
namics in dense suspensions is dominated by short-range in- 
teractions [30] . We complement our experiments and particle- 
based simulation studies by identifying a minimal continuum 
model for incompressible active flow that combines elements 
from the Toner- Tu [T5l fTHl H7] and Swift-Hohenberg [45] the- 
ories. 
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Theory and Experimental Results 

Motivation for the SPR Model. Self-motile bacteria may form 
meso-scale vortex patterns if their concentration is sufficiently 
large [3 El \M CE] • At very high volume fractions (0 > 40%), 
steric repulsion and other short-range interactions (e.g., lu- 
brication forces, flagellar bundling of neighboring cells) can 
be expected to govern physical reorientation and alignment, 
whereas intrinsic Brownian motion effects 30 become less im- 
portant in this collision-dominated high-density regime |46| . 
Chemotaxis [18] can strongly affect bacterial dynamics 
in droplets or near liquid-gas interfaces, but is less relevant 
in closed chambers as considered in our experiments. Re- 
cent direct measurements of individual Escherichia coli flow 
fields [30 suggest that hydrodynamic far-field interactions are 
negligible for bacterial reorientation, especially when bacteria 
swim close to a no-slip surface. Earlier experiments [8l l24ll25] 
on 2D films and 3D bulk suspensions also show that the aver- 
age swimming speeds of individual bacteria (typically of the 
order of 10 /xm/s in isolation 8 , 30 ) can be enhanced up to 
five times through collective hydrodynamic near- field effects. 
In the simplest approximation, however, a sufficiently dense 
bacterial suspension can be viewed as a system of determin- 
istic, self-propelled, rod-like particles with an effective swim- 
ming speed V (for B. subtilis at <j> ~ 40% we find V ~ 30 to 
100/im/s depending on oxygen concentration and boundary 
conditions). One of our objectives is to test such a minimal 
model against experiments in the limit of highly concentrated 
suspensions and to provide systematic guidance for more ac- 
curate future models. 

Non-Equilibrium Phase Diagram of the SPR Model. To iden- 
tify generic requirements for the formation of turbulent phases 
in active systems, we performed simulations of a minimal 2D 
SPR model with periodic boundary conditions (see SI Ap- 
pendix for details). In its simplest form, the model assumes 
that a rod-shaped self-propelled particle moves deterministi- 
cally in the overdamped low- Re regime with an effective swim- 
ming speed V, while interacting with the other particles by 
steric forces. Mutual repulsion is implemented by discretizing 
each rod into spherical segments and imposing a repulsive 
Yukawa force potential ~ exp(— r/A)/r, where r is the dis- 
tance, between the segments of any two rods (i.e., the decay 
length A > defines the effective diameter of a rod of length 
£). If two sufficiently long rods perform a pair-collision, this 
short-range interaction results in an effective nematic (apo- 
lar) alignment, before the rods become pushed apart by the 
repulsive force. 

Depending on the effective volume filling fraction <j) and 
the rod aspect ratio a, both defined in terms of the scale pa- 
rameter A and rod length £, the SPR model exhibits a range 
of qualitatively different dynamical phases (Fig. [I]). The nu- 
merically estimated non-equilibrium phase diagram (Fig. [l|\) 
illustrates the importance of the effective particle 'shape in 
2D: Upon increasing 0, short rods undergo a transition from 
a dilute state (D), with little or no cooperative motion, to 
a jammed state (J); this transition can be identified by the 
mean square displacement per particle, which drops off nearly 
two orders in magnitude along the transition curve. By con- 
trast, very long rods (a > 13) do not jam at moderate filling 
fractions but exhibit swarming (S) behavior and large spa- 
tiotemporal density fluctuations. Generally, the transitions 
from the dilute phase (D) to cooperative motion (regions S, 
B and T) can be characterized by the Onsager overlap den- 
sity [47]. Upon increasing cj) further, very long rods tend to 
assemble in homogeneous lanes (L), corresponding to quasi- 
smectic regions of local polar order; the swarming-to-laning 
transition is signaled by a discontinuous increase in the corre- 
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Fig. 1. (A) Schematic non-equilibrium phase diagram of the 2D SPR model 
and snapshots of six distinct phases from simulations: D-dilute state, J-jamming, S- 
swarming, B-bionematic phase, T-turbulence, L-laning (see also SI Appendix, Fig. S2 
and Movies S01-S06). Our analysis focusses on the turbulent regime T. (B) Enstrophy 
per unit area fl in units (V/X) 2 for different aspect ratios a = £/X, obtained from 
SPR simulations with N ~ 10 4 to 10 5 particles. The maxima of the enstrophy 
indicate the optimal filling fraction for active turbulence and mixing at a given value 
of the aspect ratio a. Note that values cf) > 1 are possible due to the softness of the 
repulsive force (see SI Appendix for simulation parameters). 

lation length of the two-particle velocity correlation function. 
The swarming (S) and laning (L) phases adjoin a so-called ac- 
tive bionematic [18] phase (B), where vortices and extended 
jet-like structures coexist 48 , 25 ; this phase is characterized 
by large fluctuations of the local vortex density. Most impor- 
tantly for the present study, however, the SPR model predicts 
homogeneous turbulent states (T) at high filling fractions and 
intermediate aspect ratios 3 < a < 13, a range that covers typ- 
ical bacterial values (e.g., 2 < a < 4 for E. coli and 2 < a < 10 
for B. subtilis (SI Appendix, Fig. S7). The transition between 
bionematic and turbulent phase is also signaled by the veloc- 
ity distribution, correlation functions and density fluctuations 
(SI Appendix, Fig. S3, S4). 

Homogeneous Turbulent Phase in the SPR Model. A typical 
turbulent flow state as found in the simulations, and the as- 
sociated (pseudo-scalar) 2D vorticity field uj — d x v y — d y v x , 
are shown in Fig. [5] The mean local flow field v(t,r) at 
time t and position r was constructed by binning and aver- 
aging individual particle velocities, using a spatial resolution 
similar to that in our experiments (SI Appendix). To charac- 
terize the emergence of homogeneous turbulence in the SPR 
model in terms of particle geometry a and effective volume 
fraction 0, we quantify the vortical energy through the en- 
strophy [35] [34] [27] per unit area, Q, = \ (\uj(t, r)| 2 ), where 
brackets ( • ) indicate spatial averages and overbars denote 
time averages. For slender rods (a > 3) the mean enstrophy 
Q exhibits a maximum when plotted versus the volume frac- 
tion cj) (Fig. [l£).This maximum coincides approximately with 
the transition from the bionematic to the turbulent phase; in 
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Fig. 2. Experimental snapshot (A) of a highly concentrated, homogeneous quasi- 
2D bacterial suspension (see also Movie S07 and Fig. S8). Flow streamlines v(t, r) 
and vorticity fields cu(t,r) in the turbulent regime, as obtained from (B) quasi- 
2D bacteria experiments, (C) simulations of the deterministic SPR model (a = 5, 
(j) = 0.84), and (D) continuum theory The range of the simulation data in (D) 
was adapted to the experimental field of view (217 /im X 217 /J,m) by matching the 
typical vortex size (scale bars 50/im). Simulation parameters are summarized in the 
SI Appendix. 

a bacterial suspension, it corresponds to the optimal concen- 
tration for fluid mixing. Typical aspect ratios of bacterial cell 
bodies in our experiments lie in the range 2 < a < 10 (mean 
6.3 =b 1.2; see SI Appendix, Fig. S7). Hence, homogeneous 
bacterial turbulence should be observable in 2D for a broad 
range of filling fractions. 

Experiments. We test the T-phase of the SPR model against 
experimental observations of B. subtilis at high filling frac- 
tions (0 > 50%, see Materials and Methods). In contrast to 
recent investigations of bacterial dynamics in 2D free-standing 
films [8], on 2D surfaces [44, 49 and open 3D bulk suspen- 
sions [7, 24 , 25 , 18 , bacteria were confined in closed microflu- 
idic chambers to minimize oxygen gradients that may cause 
anisotropic streaming of the oxytactic B. subtilis bacteria [2] . 
To study the effects of dimensionality and boundary condi- 
tions, experiments were performed with two different set-ups: 
quasi-2D microfluidic chambers with a vertical height H less 
or equal to the individual body length of B. subtilis (~ 5 /im), 
and 3D chambers with H « 80 /im (Fig. S6, S8 and Movies 
S07-S10). To focus on the collective dynamics of the microor- 
ganisms rather than the solvent flow [241 150] , we determined 
the mean local motion of B. subtilis directly using particle 
imaging velocimetry (PIV; see also SI Appendix). A typical 
snapshot from a quasi-2D experiment is shown in Fig. [2]\. As 
evident from the inset, local density fluctuations which are 
important in the swarming/flocking regime [511 149] become 
suppressed at very high filling fractions (SI Appendix, Fig. 
S5). The corresponding flow fields (Fig. [2p and SI Appendix, 
Fig. S8) were used for the statistical analysis presented below. 

Continuum theory. The analytical understanding of turbu- 
lence phenomena hinges on the availability of simple, yet suffi- 
ciently accurate continuum models [27] . Considerable efforts 
have been made to construct effective field theories for ac- 



tive systems [15] \W\ \T7\ \M EH E2 \M [53j [54] but most of 
them have yet to be tested quantitatively against experiments. 
Many continuum models distinguish solvent velocity, bacterial 
velocity and/or orientational order parameter fields, resulting 
in a prohibitively large number of phenomenological parame- 
ters and making comparison with experiments very difficult. 
Aiming to identify a minimal hydrodynamic model of self- 
sustained meso-scale turbulence, we study a simplified con- 
tinuum theory for incompressible active fluids, by focussing 
solely on the experimentally accessible velocity field v(t,r). 
By construction, the theory will not be applicable to regimes 
where density fluctuations are large (e.g., swarming or flock- 
ing), but it can provide a useful basis for quantitative com- 
parisons with particle simulations and experiments at high 
concentrations. 

We next summarize the model equations; a detailed moti- 
vation is given in the SI Appendix. Since our experiments 
suggest that density fluctuations are negligible (Fig. 2A) 
we postulate incompressibility, V • v = 0. The dynam- 
ics of v is governed by an incompressible Toner- Tu equa- 
tion [15] [161 [T7] , supplemented with a Swift-Hohenberg-type 
fourth-order term [45] . 

(d t + Xov • V)v = -Vp + Ai Vv 2 - (a + f3\v\ 2 )v + 

r v 2 v-r 2 (v 2 ) 2 v, [l] 

where p denotes pressure, and general hydrodynamic con- 
siderations [52] suggest that Ao > l,Ai > for pusher- 
swimmers like B. subtilis (see SI Appendix). The (a,/3)-terms 
in Eq. Jl| correspond to a quartic Landau-type velocity po- 
tential p3[l6][T7]. For a > and ft > 0, the fluid is damped 
to a globally disordered state with v = 0, whereas for a < a 
global polar ordering is induced. However, such global polar 
ordering is not observed in suspensions of swimming bacteria, 
suggesting that other instability mechanisms prevail [53]. A 
detailed stability analysis (SI Appendix) of Eq. implies 
that the Swift-Hohenberg-type (ITo, r2)-terms provide the sim- 
plest generic description of self-sustained meso-scale turbu- 
lence in incompressible active flow: For To < and T2 > 0, 
the model exhibits a range of unstable modes, resulting in 
turbulent states as shown in Fig. [2p. Intuitively, the (To, Vz)- 
terms describe intermediate-range interactions, and their role 
in Fourier-space is similar to that of the Landau-potential 
in velocity space (SI Appendix). We therefore expect that 
Eq. describes a wide class of quasi-incompressible active 
fluids. To compare the continuum model with experiments 
and SPR simulations, we next study traditional turbulence 
measures. 

Velocity Structure Functions. Building on Kolmogorov's sem- 
inal work [55], a large part of the classical turbulence 
literature [27J El EH [37J El dOj SJ focuses on iden- 
tifying the distribution of the flow velocity increments 
Sv(tj r, R) = v(t, r + R) — v(t, r). Their statistics is com- 
monly characterized in terms of the longitudinal and trans- 
verse projections, 8v\\ = R- Sv and Sv± = T • Sv, where T = 
(eijRj) denotes a unit vector perpendicular to the unit shift 
vector R = R/\R\. The separation-dependent statistical mo- 
ments of Sv\\ and 5v± define the longitudinal and transverse 
velocity structure functions 

Sft i± (R) := ((8v n , ± ) a ) , n = l,2,.... [2] 

These functions have been intensely studied in turbulent high- 
Re fluids 27, 34 , 35 , 41 , but are unknown for active flow. For 
isotro pic steady-state turbulence, spatial averages ( • ) as in 
Eq. PJ become time-independent, and the moments Sn }± 
reduce to functions of the distance R — \R\. 
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Fig. 3. Velocity statistics of self-sustained turbulent phases in active suspensions. (A) The marginal distributions of the normalized Cartesian velocity components 
[vi — (vi)]/[(Vi) — { v i) 2 \\= x y are approximately Gaussian (thin grey line) for experiments, SPR model and continuum theory. (B) The distributions of the longitudinal 
and transverse velocity increments Sv \\,±, normalized by their first and second moments are shown for three different separations R. (C) Longitudinal and transverse 

velocity structure functions Syf j_ normalized by (v 2 ) ^ . The maxima of the even transverse structure functions S 2 ^ reflect the typical vortex size R v which is significantly 
larger in the 3D experiments. Experimental and theoretical data points are spatio-temporal averages over two orthogonal directions in (A) and (B), and all directions in (C), 
yielding a typical sample size > 10 6 per plotted data point in (C). Histograms and structure functions for quasi-2D (3D) curves were obtained by combining PIV data from 
two (fifteen) movies, respectively, representing an average over 2x1000 (15x300) frames. Simulation parameters are identical to those in Fig. [2] and summarized in the SI 
Appendix. Errorbars are smaller than symbols. 



Velocity distributions, increment distributions and struc- 
ture functions for our numerical and experimental data are 
summarized in Fig. [3] For the SPR model, the velocity statis- 
tics can be calculated either from the raw particle data or 
from pre-binned flow field data. The two methods produce 
similar results, and Fig. [3] shows averages based on individual 
particle velocities. Generally, we find that both the 2D SPR 
model and the 2D continuum simulations are capable of repro- 
ducingthe experimentally measured quasi-2D flow histograms 
(Fig. K3K,B) and structure functions (Fig. p3C). The maxima 
of the even transverse structure S^ 1 signal a typical vortex 
size R v , which is substantially larger in 3D bulk flow than in 
quasi-2D bacterial flow. Unlike their counterparts in high-Re 
Navier-Stokes flow [271 134], the structure functions of active 
turbulence exhibit only a small region of power lawgrowth 
for £ < R <C R v and flatten at larger distances (Fig. [§□). 

Velocity Correlations and Flow Spectra. The energy spectrum 
E(k), formally defined by (v 2 ) — 2 J °° E(k)dk : reflects the ac- 
cumulation of kinetic energy over different length scales. By 
virtue of the Wiener-Khinchine theorem [27], E(k) can be esti- 
mated by Fourier-transformation of the equal-time two-point 
velocity correlation function, yielding in d dimensions 

E d {k) = J d d R e~ ihR {v(t, r) ■ v(t, r + R)), [3] 

where C2 = 2tv and C3 = 47T. Normalized velocity correla- 
tion functions (v(t, r) • v(t, r + R)) and spectra Ed(k) for our 
data are summarized in Fig. [4] The crossover from positive 
to negative correlations indicates again the typical vortex size 
R v , in agreement with Fig. [3p and previous findings for open 
3D bulk systems [TlfTH]. 

In bacterial suspensions, the microorganisms inject ki- 
netic energy on small scales R ~ £, setting the upper bound 



kg = 2n/£ for the spectral range of the bacterial fluid. For 
both experiments and simulations, we observe turbulent vor- 
tices on scales R > £, which formally correspond to the energy- 
inertial range k < kg in classical 2D turbulence 34, 35 . Our 
experimental and numerical data suggest asymptotic power 
law scaling regimes for small and large k- values (see Fig. [4^3) , 
but the power-law exponents differ from the characteristic 
/c" 5/3 -decay of 2D Kolmo gorov-Kraichnan turbulence [39] : 
see discussion below. The spectra for the 2D continuum 
model and the quasi-2D bacteria experiments are in good 
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Fig. 4. Equal-time velocity correlation functions (VCFs), normalized to unity at 
R = £, and flow spectra for the 2D SPR model (a = 5, 4> = 0.84), B. subtilis 
experiments, and 2D continuum theory based on the same data as in Fig. [3] (A) 
The minima of the VCFs reflect the characteristic vortex size R v [48]. Data points 
present averages over all directions and time steps to maximize sample size. (B) For 
bulk turbulence (red squares) the 3D spectrum Es(k) is plotted (kg = 2tt/£), the 
other curves show 2D spectra E^ik)- Spectra for the 2D continuum theory and 
quasi-2D experimental data are in good agreement; those of the 2D SPR model and 
the 3D bacterial data show similar asymptotic scaling but exhibit an intermediate 
plateau region (spectra multiplied by constants for better visibility and comparison). 
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agreement, both showing large- k scaling with approximately 
E(k) - k~ 8/s and small-*; scaling with roughly E(k) ~ /c +5/3 . 
The asymptotic spectra for the 2D SPR model and the 3D ex- 
perimental data look qualitatively similar, but do also exhibit 
an intermediate plateau region which indicates that kinetic 
energy is more evenly distributed over a range of scales. 



Discussion and Conclusions 

SPR Model vs. Experiment. The deterministic SPR model 
provides a simplified description of the bacterial dynamics, 
as it neglects not only elastic properties of flagella and cell 
body, but also hydrodynamic interactions and orientational 
fluctuations due to intrinsic swimming variability and ther- 
mal effects [301146] . Notwithstanding, at high concentrations, 
such a minimal model reproduces remarkably well the flow 
velocity distributions and the structure functions from our 
quasi-2D B. subtilis experiments and the 2D continuum sim- 
ulations (Fig. [3]). This implies that hydrodynamic interac- 
tions per se are not required for the formation of self-sustained 
turbulence in dense suspensions - self-propulsion, a rod-like 
shape and volume exclusion interactions are sufficient (this 
raises the question whether the optimization of collective be- 
havior may have been a factor in the evolution of bacterial 
shapes). However, to achieve a better quantitative agreement, 
particle-based future studies should focus on more realistic 
models that account for hydrodynamic near-field interactions 
and intrinsic randomness in bacterial swimming [30] . The ex- 
perimental results presented above provide a benchmark for 
evaluating such microscopic models [56] . 

Continuum Model and 'Universality'. The good agreement of 
the structure functions, spatial and temporal flow correlations 
(see also Fig. S9), and spectra from the 2D continuum the- 
ory with those from the quasi-2D experiments suggests that 
this theory could be a viable model for dense suspensions. 
Since the instability mechanism in the continuum theory arises 
from a generic small- wave number expansion in Fourier-space 
(see SI Appendix) , that is analogous to the Landau-expansion 
in order-parameter space for second-order phase- transit ions, 
we expect that the model applies to a wide range of quasi- 
incompressible active fluids. This would imply that meso- 
scale turbulent structures in these systems share 'universal' 
long- wave length characteristics. We note that the theory as 
formulated in Eq. Q only accounts for leading terms up to 
fourth-order and, therefore, becomes inaccurate for large ve- 
locities and wave numbers (see tails in Figs. 3A,B and 4B). 
Nevertheless, this continuum model appears to capture the 
main statistical and dynamical features of the experimental 
data. Important future challenges include the analytical pre- 
diction of active flow spectra from Eq. detailed numeri- 
cal studies of 3D bacterial bulk flows, andcomparisons of our 
experimental and numerical data with Q-tensor models and 
other multi-order parameter theories [2] [TBI El EB E2] • 

Dimensionality, Boundaries and Hydrodynamic Interactions. 

The quasi- 2D experiments allow us to compare with 2D sim- 
ulations that come close to experimental system sizes. Free- 
standing thin films 8 and bacterial mono-layers on open sur- 
faces [441 149] , which may be more prone to intrinsic insta- 
bilities and external fluctuations, provide an alternative, but 
non-equivalent realization of a 2D bacterial fluid. The cru- 
cial difference between freestanding 2D films and our closed 
quasi-2D set-up is that the presence of no-slip boundaries in 
our experiments suppresses hydrodynamic long-range inter- 
actions between bacteria due to cancellation effects from the 



hydrodynamic images: An isolated dipole-like swimmer (as E. 
coli [30] and, most likely, B. subtilis) creates a stroke-averaged 
far-field flow that decays as ~ 1/r 2 with distance r in a 3D 
fluid. When the same swimmer moves parallel to a nearby 
solid surface in an otherwise semi- infinite fluid, the flow com- 
ponents parallel to the boundary decay faster ~ 1/r 4 [30] . If, 
however, the swimmer is closely confined between two parallel 
no-slip walls, as in our quasi-2D experiments with H ~ 4 /im, 
then the flow field becomes exponentially damped at distances 
\r\ ^> H [57] . By contrast, in free-standing 2D films the 
flow field generated by an isolated microorganism has a much 
longer range ~ 1/r 58, 22 , suggesting that hydrodynamic 
interactions could play a more important role for collective 
behavior in these systems [8] . The fact that the typical vortex 
size in 3D is larger than in quasi-2D could indicate stronger 
short-to-intermediate-distance hydrodynamic coupling in 3D 
bulk flow; it would therefore be interesting to perform a sim- 
ilar analysis for thin-film data [8]. Generally, however, we 
expect hydrodynamic far-field interactions to be less impor- 
tant for the dynamics in very dense suspensions due to mutual 
hydrodynamic screening [59 and the small magnitude of bac- 
terial flows fields [30 , but they could act as a destabilizing 
noise [541 160] . 

Low- Re vs. High- Re Turbulence. Conventional high- Re tur- 
bulence arises from energy input on large scales (e.g., stir- 
ring or shearing). In 3D flow the injected energy is redis- 
tributed to smaller scales via an energy-inertial downward 
cascade with E 3 ~ k~ 5/s [27]. In 2D films, due to the 
suppression of vortex stretching [35] [34] , there can be both 
an energy-inertial upward cascade with E2 ~ k~ 5 ^ 3 and an 
enstrophy-transfer downward cascade with E2 ~ k~ 3 [39] . 
Remarkably, viscoelastic polymer solutions can exhibit tur- 
bulent features (e.g., spectral power law scaling) at Reynolds 
numbers as low as 10 -3 , facilitated by a slow nonlinear re- 
sponse to external shear due to long intrinsic relaxation times 
of the polymers [61] [62]. Our simulations and experiments 
suggest asymptotic spectral power law decays towards the 
bacterial energy injection scale kg — 2tv/£ that resemble the 
energy-inertial regime of classical turbulence but, due to vis- 
cous damping by the low- Re solvent, extend over a smaller 
range of length scales (roughly up to 10£). The latter fact is 
reminiscent of viscoelastic turbulence [61], although the un- 
derlying physical mechanisms are very different. 

In conclusion, bacterial or, more generally, self-sustained 
'active turbulence', shares some qualitative characteristics 
with classical turbulence on small scales while differing on 
larger scales. Our detailed statistical analysis shows that, as 
with inertial turbulence, a complete quantitative understand- 
ing of turbulent behavior in active systems poses a challenging 
task. The combined experimental, theoretical and numerical 
results presented here may provide both qualitative and quan- 
titative guidance for future studies that aim at identifying the 
basic principles of dynamical self-organization in living fluids. 
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Materials and Methods 

B. subtilis cells (wild type strain 168) were streaked from a -80° C stock onto an 
LB medium plate containing 1.5% agar. The plates were incubated at 37° C for 12 
h. A single colony from the plates was used to inoculate an overnight culture in 
Terrific Broth (Sigma), which was then back-diluted 1:200 into 50 mL of fresh tryp- 
tone broth, and grown at 37° C on a shaker to mid-log phase. The culture was 
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then concentrated 400 X by centrifugation at 4000 Xg for 3 min, and the pellet was 
resuspended by gentle vortexing, to not shear off the flagella. The concentrated cul- 
ture was loaded into a polydimethylsiloxane (PDMS) microfluidic device, which was 
then sealed to reduce background fluid motion. The microfluidic device consisted of 
cylindrical measurement chambers (radius 100 /im, height 4 /im for quasi-2D mea- 
surements, and radius 750 /im, height 80 flm for 3D measurements). The samples 
were imaged in bright field with a 40X/NA 1.4 oil immersion objective on a Nikon 
Tl-E microscope. Images were acquired at 40 fps in 2D (camera: Pike, Allied Vi- 
sion Technologies), and 100 fps and 200 fps in 3D (camera: Phantom v9.1, Vision 



Research). Compared with measurements in quasi-2D chambers at the same frame 
rate, the vertical superposition of bacteria leads to a reduced image quality in 3D 
samples; we therefore recorded the flow in 3D suspensions at a higher frame rate. For 
the 3D measurements, we imaged at the bottom and in the middle of the chamber, 
while for the quasi-2D measurements, we imaged in the middle of the chamber. A 
detailed description of the theoretical models and numerical methods is given in the 
SI Appendix. Raw data and additional experimental movies can be downloaded from: 
http://damtp.cam.ac.uk/user/gold/datarequests.html 
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